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ABSTRACT 


A statistical approach to sound propagation is considered 
in situations where, due to the presence of large gradients of 
properties of the medium, the classical (deterministic) treatment 
of wave motion is inadequate. Mathematical methods for wave 
motions not restricted to small wavelengths (analogous to known 
methods of quantum mechanics) are used to fcmulate u wave theory 
of sound in nonuniform flows. Nonlinear transport equations for 
field probabilities are derived for the limiting case of noninter- 
acting sound waves and it is postulated that such transport 
equations, appropriatly generalized, mav be used to predict the 
statistical betvwlor of sound in arbitrary flows. 
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T. imHouucTTOtt 


In this work we address ourselves to the nonelasslcal treat* 
went of sound as a wave phenomenon. By "nonelasslcal treatment" 
we shall mean an analysis based on mathematical methods developed 
for the description of quantum systems in cases where wave proper- 
ties change appreciably over one wavelength. The objectives of 
this paper are to formulate a quantum-like stochastic theory of 
sound of arbitrary intensity propagating through a perfect gas 
in an arbitrary irrotational flow. 


A classical approach to wave analysis is valid only If 
|VX|«l, where ' * wavelength. Acoustics of inhomog,eneous moving 
fluids concerns itself usually with the zero-wavelength approxi- 
mation, V? » 0. An extension of the classical approach in the 
spirit of the WKB (Wentzel-Kramers-Brillouin) approximation, 

Kentzer (1975a), resulted in a dispersion relation valid for small 
but not necessarily negligible wavelength. However, situations 
where jVX| - 0(1) are not uncommon in high speed flows. This 
may be demonstrated on the basis of a plane wave propagating through 
a fluid which moves with velocity U. Tf the speed of sound is a, 
then the frequency is w * (U+a)k and we have along the wave tra- 
jectory (here k * 2x7 X » wavenumber) 


dk 

dt 





dx 

dt “ 3k 


(lf+a) 


so that. In this case, 


dX 

dx 



Thus the gradient of wavelength is equal to the fractional change 
in the propagation velocity, (U+a>, occurring over one wavelength 
Therefore a cieualcal approach is valid if \<<|(U+a)/V(U+a)j . 


There are manv applications, especially in the field of aerodynamic 
noise, where this inequality is violated. For example, when a sound 
wave props pates upstream a small change in (U-a) mav produce a large 
fractional change in (U-a) especially when U 'v a or at near-aonic 
conditions. Another situation demanding nonclasaical treatment is 
one in which sound propagates across boundary layers or shear layers 
fOv then 3U/3x is large. A very special situation arises in turbu- 
lent flows where, due to the "graininess" of the medium, 3U/3x is 
very large almost everywhere and may change sign several times in 
the distance of one wavelength. When vorticity or entropy gradients 
are present, one must consider a more general theory which would 
include the interactions of acoustic waves with vorticity and entro- 
py waves. Such a theory, based on the full Navier-Stokes equations 
of fluid mechanics, was developed by Kentzer (1975b). 

In this paper we limit ourselves to sound propagating through 
arbitrary potential (irrotational) flows of a perfect gas. The 
organization of the paper is as follows. First, we determine the 
normal modes of the linear part of the equation for the fluctuations 
around the instantaneous local mean flow, Sec, II. The normal modes 
which form a doubly infinite set of orthogonal basis vectors, are 
used to form an integral representation of the solution of the non- 
linear equation. The exact nonlinear equation is satisfied Only 
locally on the average under the assumption that the phase of the 
fluctuations is a random function. Such a function is then modelled 
so as to preserve a wave-particle duality. A complex characteristic 
function t|/ is introduced in Sec. Ill, and an equation in the form 
of a nonlinear, convective and dissipative Schrodinger equation for 
the characteristic functions is derived treating local mean flow 
conditions as constant. The Schrodinger-type equation is then 
recast In Sec. IV hack into a hydrodynamic form by Madelung's 
transformation In which the new dependent variables are now the 
absolute value of the magnitude and the phase of the characteris- 



tic functions. The square of the magnitude, l*I'|% is then inter- 
preted as a probability density and the gradient of the phase, when 
multiplied by a constant with units of action density, is inter- 
preted as a probability diffusion velocity. In Sec. V we show 
how ensemble averages (moments) of products of wavenumber components 
may be defined using tbe characteristic functions, and how the 
moments determine the wavenumber distribution functions. The dis- 
tribution functions, in turn, are used to calculate local ensemble 
averages including the nonlinear wave interaction terms. Examples 
of averages of quantities of interest are given. Conclusions are 
summarized in Dec. VT. Tinallv, in tbe Appendix, we factor out 
the linearized convective wave equation for the two orthogonal 
eigensolutions into two linearly independent first order equations 
in the characteristic form. 


-i,- 


II . EXPANSION IN TERMS OF NORMAL MODES 
OF THE LINEARIZED WAVE EQUATION 

cole flow of a baro tropic fluid governed 

0 

- 0 , 

2 

p(p) « mass density, a*'" Yp/o, Y “ ratio 

of specific heats. 

It will be convenient to introduce a velocity potential 

the velocity vector u is given by u * Vi*, and a pressure 
P defined by P * /dp/o(p) * a~/(y-l). Thus 

u*Vu + VP « 0 
u*VP + (y-l)PV*u » 0, 

~ + 74>*W<1> + Vp - 0 

dt 

*ip 

9t + + (Y-l)PV*V-i - 0. 

We introduce now the mean (ensemble average) and fluctuation 
values, 4>» $ + <&, P » P + p*. The result is 

+ V$*VV? +■ VP] + + V$*Wd> 4- + VP'] 

dt a Z 

+ tvd>*vvd>] « o 

t~+ V$*VP + (y-l)PV* V$] + [||V V$-VP'+ V$*VP 


such that 
potential 

+ 

It 

3P . 


or, 


?t 


av* 


Consider an iscntr 
bv the equations 

+ i Vp » 
d 


^ i- u*Vp + oa^V«u 


where p * pressure, p » 


+ (Y-1)(PV*V4> + P'V-V?)] + [V4>- VF'-*- (y-1)P*7«V*] - 0 


1 


ii 


* 
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By definition, the nvernp.es of the fluctuations vanish. We may, 
therefore, average the above count Iona to obtain differential 
equations for the mean motion: 

fflft + y$.v$5 + VP + (V‘**Wf) » n Cl) 

3t 

if. + v5*yp + (y- 1 ) P ! 7 • 7? + r.'i •'.'P') + fv-DCP'^’V*) - 0 (2) 


Subtracting Eqs. (1) and (2) we obtain the differential equations 
for the fluctuations: 


+ u*V*?A + VP' - [ • 77* ] + { ( r M- ♦ V 1 '* ) - (7»b • W<t> ) ] 
Dt 

- h + »i 


U*VP'+ (y-DPV*?-!. «= [ • VP- (y-1 ) P ’ V • u l+ 1 (V4> *7P ’ ) 
3t 

- . !,-p * + (y-l) { (P '?.~Vrf) -p ' V • Vrf } ] 

- I, + K 


O) 


(A) 


where U « - mean velocity vector. In Eqs. (1) and (2) the 

averages of the products of the fluctuations are in the nature of 

Reynolds stresses. If the energy of the acoustic fluctuations is 

small, one may neglect the quadratic terms in (1) and (2) to 

facilitate computations, for then Eqs. (1) and (2) uncouple from 

Eqs. (3) and (A), and the mean flow becomes independent of the 

fluctuations. In Eqs. (3) and (A) we have separated by square 

brackets terms linear in the fluctuations, 1^, and terms quadratic 

in the fluctuations, B „ These terms will be referred to as the 

wave interaction terms. We note he^e that I vanish if the mean 

flow is uniform, and B vanish if the mean flow is uniform and 

a 

the squares of the fluctuations are negligible, a ** 1,2. 

We digress now and consider only the linear parts of Eqs. (3) 

and (A), that Is, we set I •* B - 0, U - constant, P - constant. 

a ct 



1 


II 


I 


a 


»h- 


In the linear case we mav take as the basic solution 

. i (x*k“ut) ,, i « x»k"j>t 1 

* - t*e , I’ » t' c i; , 

substitute in Fqs. (31 and (A), take a scalar product of bq. (3) 
with 0 and solve for the amplitudes ■!*« and r„: 


| u i tel i' u k V: k , i , ‘ ih i u i 


... 




-(Y-l}Pk k k k ,-l«M -iy: k > 


Po. 


tSJ 


Nontrivial solutions of Eq. (5) exist if and onlv if the determinant 
of the coefficient matrix vanishes. Setting the determinant equal 
to zero given the familiar dispersion relations for sound, either 
in the quadratic or linear form: 


- ? 


_ > i a 


or. 


(to ~ U*k) » (Y-llPk" <* a w k“ 


u • U*k + ak « U*k + c ak, c, “ 1, c. * -1 . 

— >t 1 2 


(6) 

(7) 


In the above, k * wavenumber vector, k * |k|. To the two distinct 
roots there correspond two eigenvectors, 


1J 


{1,-11. R-, « U, IK 


•1 




uj Pi “ 'nP 


f 1 a - R 
\0 a 4 R * 


The amplitudes <b 0 and P„ for the two modes tt « 1,2 are then given by 


<t> 0 - ~r H , , P„ - -iA H 
ak al* a a2 


Superposing the solutions and assuming a continuous spectrum, 
we may write 

A (k) is^(k) is n (k) 

$ -if — R al e - dk, ?' - -lE/A^CkJR^e dk 


ak 


0t2 


( 6 ) 


ot-2 


-^Z dk i dk 2 dk 3’ ^ ^ ’ s a “ ** k ” W a t ’ 

ec**l 


where . dk 




f 
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Returning now to the nonlinear equations for the fluctuations, 
Eqs, (3) and (4), we mav expand tv* dependent variables using the 
doubly infinite set of orthogonal vectors R^e ns hasia vectors 
and use the expressions (6) provided that the wave antplitudes A^ 
are allowed to depend on space and time, A f * A (x,t,k). Equations 
(3) and (4), multiplied by various powers of A^Tt^e and integra- 
ted over k, would generate an infinite set of integro-differential 
equations for amplitude correlations of various orders. Due to 
the nonlinear terms B in (31 and (4), the infinite set would be 
open, i.e., there would be always more unknown correlations appear- 
ing in the interaction terms than the number of equations. In 
order to avoid generating an open hierarchy of correlation equations 

we assume that the amplitudes A vnrv siowlv in space and have a 

o 

random phase, We will set 

A(i t t,4«A(k)e’ l '3 t (7) 

Cl Ct 

where is a complex random function. The assumption (7) amounts 
to modifying the frequency w by an addition of the random function 
With this modification of frequency, substitution of expressions 
(6) Into Eqs. (3) and (4) gives 

S/ l^a A a R al exptl(5i,C ~ V ~ V )]dR " f l + ®l 

- £/ to A R _exp(i(x‘it - to t - uj t)]dlc - I_ + B 0 . 
a a a/ ot a i l 

We now multiply the first of the above equations by 

aj~,AgRg^exp [-i(xMt-togt-togt) 3 , and the second by Agl^exp [-iCx't* 

-(Dgt-53gt) ] , subtract, and integrate with respect to k’ . Due to the 
orthogonality of the eigenvectors we have 

2/Wg!Ag| 2 dk - a/ (fj+55^) •^, , Ag(k' jRg^xp [-iCx^k'-Wgt-Wgt) ]dk' 

- (*2 + ®2^ A B (k')Rg2 ex P ["1 Jdfc' 


( 8 ) 


where A* ■ complex conjugate of A (V , oi* » w(k'), to' " w(k'). If the 
amplitudes A^(k) are known, even approximately, as functions of k, 
the Indicated integrals mav he evaluated subject to appropriate 
resonance conditions. 

We remark here thut the approximation (7) accomplishes three 

things at the same time, l a , we satis fv formallv the fluid dynamical 

equations for the fluctuations, Fqs, (3)-(4), by considering them 

to be defining equations for £o , 2 a , we can account for the nonlinear 

wave interactions in the average at a given point and instant of 

time, and, 3", we avoid the closure problem by not using Eqs.(3)-(4) 

as kinematical equations fo? v bt time evolutiun of the amplitudes 

A and their correlations, 
a 

Since I and B arc, respectively, linear and bilinear in the 

a a 

fluctuations <t> and P’, it is evident that ib^ will have separate con- 
tributions from I (two-wave interactions) and from B (three-wave 

fX 

interactions). Due to the products of exponentials appearing under 
the Integral signs in Eq.(8) the wave interactions ore interpreted 
as two— and three-wave resonances. The resonance conditions, ana- 
logous to equations of conservation of quasi-particle momentum and 
energy, are of the form 

ft ■ It’, Iji (?) + W (It) * Wrt(It’) + _ C K. * ) 

a a d o 

for pairs of interacting waves, and 

it - k'+it", w^(fc)+to a (k) » t a 8 (k , )^(k’)4 <l)Y (k 1, )^ Y (k’') 
for triads of interacting waves, nt,R,Y " 1,2. 

A preliminary studv of wave interaction integrals in Eq. (8) 

Indicates that, in absence of randomness (w^ * 0), only waves with 

parallel wavenumber vectors may interact and that an addition of 

even the smallest viscous terir.B would eliminate all interactions 

among acoustic waves of finite wavelength. However, the randomness 

of the medium (to ¥ 0) makes the resonance an effective mechanism 
a 



h 


Ji 


M - 


for the redistribution of the spectral energies. 


In order to solve to* (8) for 6,-, metre easily we will assume 
the wavenumber dependence of fly,. T* 16 lowest order nontrivial poly- 
nomial approximation to <i> (it) Is a quadratic expression, 


w^(k) 


<: m - T> 

r .2. 


% - t‘V k * 


(9) 


which expression has an advantage of being analogous to a Hamiltonian 
quadratic In momenta, at*d the advantage of affording, an interpreta- 
tion of wave resonances as interactions among quasi-particles (wave 
packets) . 


With the expression (8) for S.,fk) we reduce the integral Eq.{8) 

b 

to an algebraic equation for the complex function (r,^ - if )i 


-if . 

D 


(a/ fTj+Bj) *|.A*(k' iK^expt-i (x4-t^t-o.d t) Jdk* 

- /a 2 +H 2 )Ajd?‘ >l' K ,expf-i Cx*k-wj’t-S^t) ]dk* > 
*[2/k 2 f i; ,(k)dkf 1 

l ! 1 


where fg(k) r 



wavenumber distribution function of the 6- mode. 



II 




- 1 1 ‘ » - 


m. atAUAmwtsm: ri:‘:rn-i‘is 


Is. 

Let tl» B /A dxp [ i t-ru tlldk * /A e l dk 
a it n « n 

where a * x*C»w t-to t is the phase. If A Is a function of k only, 

>x a a i 

then 

M’ - ^ 9 ' t - 

~-J* » -i / (>.:> +fj* I A £k ) e Aik 
(it an ft 

Is , is 

*r»l> * ij’kA e '*JV; , i r *I» ® - /k A e 'dk, 
a a m u 

Consequently, to the dispersion relation 
[wd-6h*(c-i£.)k fc ~n«Fp - • *> 


there i;** responds the differential equation 


U(|r + 5*V) + (r.-iOvV* + iC't* 

o t 


which has to be satisfied by the characteristic function <k for 
both modes rt * 1 and it * 2 , 


The obiections to the use ot the wave equation for •[» are the 
J u 

need for two initial conditions for the charaeterisitc function of 
each mode and the difficulty of separating the magnitude and phase 
of the chorcteristic function in an equation of second order. When 
£3 * 0, the equation reduces to the wave equation. As shown in the 
Appendix, if U and a are kept constant the wave equation may be 
replaced by a svstem of two equations of first order in a characte- 
ristic form: 

[“+ (U + c^an) * 0, Cj - 1, c 2 - -1, 

with the arbitrary unit vector n being the space component of the 
characterisitc normal. This result suggests the following heuris- 
tic approach. We first factor out the dispersion relation into 


a 


ft 
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two factors, m * U»l + c ak *■ ('. ~t' f )k\ and write k ** n*k where 
a a >i *t 

n “ £/k * unit vector in the direct ton of k. Thu*, the differential 
operator corresponding to k is -in*”. In anticipation of a genera- 
lization of t tie results to uonuniforr? flows, we will first symmet- 
rize the dispersion relation. Subject to verification by experiment, 
we shall take 

w. - ~(U*k + k*P) * uti *v + n*ka) - fk *(C -if, )k] “ 0 

U w- *• (I tt 


to which there corresponds the differential equation For the charac- 
teristic function of the w-mode: 


*► 5 tU^ V 1 k -c ( [an *'?*_■ i + n # V(a ; l' ) ] 

4 ■ ♦ i (. -it. r-’-p i - ». 

ct Ci u 


This may be rewritten in the form of a nonlinear, convective and 
dissipative Schrodinper equation, 


3ti* 

7“ a + m 4 c a*") * V7.f,t g r{ , ) * 17. (r 0(1) ) 

et a a it u a a 

- ^ I' k’ • [u 4 c an]. (10) 

2 Cl Ct 


We observe that C plavs the role of a diffusion coefficient while 
a 

C is an analog of the Planck constant ft * h/2tr. The units of C 
are (frequency) * (length.) ** (energy) *(time)/ (mass) , or ” specific 
action (action per unit mass). By analogy to ii, £ will determine 
the local proportionality (scale) factor of wave and mechanical 
attributes in the wave-particle duality. Classical limit is ob- 
tained for and E approaching zero. This limit corresponds to 
I ■ - 0, the case of steadv uniform mean flow with negligibly 

small amplitudes of the fluctuations. 


TV. PROBABILITY Tr Vr ;.•>>) r EOUATJnNS 


We shall put now the Schrodtnp.er-tvpe equation (10) for the. 
complex characteristic function # Into Madelung'e hvdrodynamicnl 
form, which form remains nonlinear even when U and a are constant. 


Let 


<!> 


P^e ' l where P 

Ot i.'t 


fji >i>* 
ft <t 




- probabllltv density, and where S ® /V»dx, V»V + 0 in general. 

Note that the units of V are (length) 1 and that * velocity. 

Thus V will be referred to as the "probability velocity" with 
5 being the scale factor. 

Substituting the polar form i * P*e “ into Eq, (in), dividing 
by i]i and separating the real and imaginary part 6 1 one obtains 
two nonlinear equations 


3P o 

r- +(W c an+2£ V)-9P, - V* (f \*P ) 

3 1 a a ot a n 

- {?• (U+c c a£)+Z e »* CC a V)+f Vj t j(t^.MP a ) W } }P a (11) 


3S 

^- a +(U+c an) ‘VS - *?*[E V + ±, : . VtnPJ + £ V ♦ VlnP 
3t a a a ct 2 a a ot a ot 


+ c, f~(?6tP ) ? - v 2 1 . 

Ot 4 rt n 


( 12 ) 


Equation (11) Indicates that the probability density P q is 
not conserved, is convected by the fluid and diffuses relative to 
it with the speed of sound modified bv the effects of wave inter- 
actions. The right-band-side of (11) represents a sum of dissi- 
pation and source terms. The two modes, ot = 1,2, couple only 
through the coefficients r and r .^. 

Since the phase S of the characteristic function has no 
physical meaning, and since only its gradient appears in the 
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transport equation for P , 
obtain 
OV 

^ ri +(U+c an) •'"V - W 

at U ' X 


+ «■ v •vHnP 

a i .i >x 


we take the gradient of 1-q, (II?) and 


(' V + k ' OtV ) + '. [i(’’5tP ) 2 -V J ] 
Vi « 2 'u a i\ 4 u a' 


} - V * P ( u+ c an! - t ,t+ c an)*(-2*V ) 
>x ft ti Ot 


- V *( t ?*(ii+c an) ) . 


(13) 


We note that the hvdrodvnamic transport equations (11) and (13) 
are nonlinear even in the case of a stationary medium, U * 0 
and a « constant, and that P. and V are strongly coupled. Due 
to the nonlinearitv of the hydredvnamical form • f Eqs. (11) and 
(13) we feel confident to postulate that these equations hold 
even In flows with arbitrary nonconstant U and a (in inhomogeneous 
media or nonuniform flows). 

As a consequence of the formulation, the number of dependent 

variables is eight , one scalar P and one vector V for each modfs 

<x <x 

a «* 1,2. However, the spatial distributions of P^ and will 
provide essential information on the spectral distribution of 
the fluctuations of the original two scalar unknowns. 


V. STATISTICAL. PROPERTIES OF FLUCTUATIONS 


If we define the averages of n function of k with respect 
to the distribution function f^Ot) as 

<F(k)> - /F(k)f u (k)dk / /f a (k)dk, 


then, by repeated differentiation we mav obtain averap.es of a 
product of the components of k: 


m i "n l 

<k \ .,k 
n no 2i 


< i 


Jv, 

5 i' 


a 


i"! m, nu, 
fJt 3x 1 . . . 3x 
n l n N 


1 

<]'* 


..K, 

,j tj i 


m. 


a .. 1 . 

>lx ... 3x 

"l n N 


where K - m.j + ...m fJ ■ the order of the product, and where n^.n^ 

indicates the order in which the space derivatives are taken. 

The averap.es may be written entirely in terms of and V , e.g.. 


'V ■ V <k iV 


.1 ■ l r „ 3fnP. ,, 3£nP, , .. . H 


If the distribution function f is also a function of x and 
t, the above definition of averages may be considered to hold 
locally with x and t as parameters. Since the characteristic 
function determines the averages of products of k of all orders, 
and the moments of the distribution determine the distribution, 
the theory is closed. That is, if tKx*,t.) and its derivatives 
evaluated at x, and t« are given, we can determine f(x»,t,,k) 
and evaluate the interaction integrals in F.q. (8). With C a *nd 
known, one may use Eqs. (11) and (13) to advance the solution 
to the time t#+dt. 

A finite number of moments may be used to determine s finite 
number of terms of an expansion of the distribution function f(fc). 
For that purpose the most convenient is the Gram-Charlier expan- 
sion using Grad polynomials (multidimensional Hermite polynomials). 
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Such polynomials re Introduced by H. Grad (1949a) and used In 
connection with multidimensional Gram-Charlier expansion to de- 
velop his 13-moment approximation to the kinetic theory of gases 
(Grad, 1999b). The convenience of the Gram-Charlier expansion 
lies in the fact that the coefficients of the series are the 
averages of the Hermite polynomials which may he written in terms 
of symmetrized averages of the products of k. Thus a truncated 
Gram-Charlier expansion will produce a distribution function 
having a finite number of prescribed moments. We may write 




where 


H 


+ (H 200^020' H/ fi02* m ilO fW 10l' H{ 01 1 ) 

+ (H 300^030^003^210^201^120^102 4H 021 + 8l2' H, Ul ) 

+ .... J, 

2 2 2 

.a. 3/2 ■ atk l +k 2 Hc 3 ) 3. , 2 .2 , .2 .-1 

g - (“) e * a - 2 ^ k l + k 2 + ** 3 *! * 

»ijk H ijk 


ljk ilJJk! (2a) 1+j+k * 


.i+j+k. 


H... (It) * (-1)* ^ k - — ™r — ~-£— * Hermite polynomial of 

g Ok^Dk^k^ order i+j+k , 


" <H^ k > “ average of Hermite polynomial (symmetrized) 
for example, H 20Q - <2o) 2 (k^ - ~] , - (2o) 2 [<k 2 > - 

H 110 " (2o) 2 k^k 2 , H 110 - (2<7) 2 [<kjk 2 > + <^>1/2. 

is 0 

Treating the Fourier components 4> * (A /ak)R ,e “and 

Is a a cu 

” ^a R a 2 e ” S8 tke 9 uant i tles t0 be averaged, and writing for 


the Induced velocity u. * V<t*. * ikA , we have 

ct Ot :t 


<<J>> « “ <P'> “ 0, 


0 , E/k“ 2 A A*f(k)dk 

«r> - «bA*> - K -iULO— 

a" >:/f a (k)dk 


<u 


<P 


' 2 > m <u f . U '*> " - 


1 Url A t A*f t Ck)dk x E/f 2 Ck)dk 


a 2 Iff (k)dk a 2 £/f (k)dk ’ 

a ( * 

^ jT 


, * E/A A*f (k)dk E/f‘(k)dk 7 - 

• 2 > - <P'P’*> - - » 2— - a 2 <u‘ 2 >. 


V f (k)dk t/f <k)dk 
a a 


If the pressure perturbations are small compared to the mean 
pressure, |p* | <<p » we may use the Binomial Theorem to obtain 

p* K p'/p so that with p « —4 ^ we have p* * pP'. Then 
r Y-l ’ Y-l ~ 


<P’ 


,2 > « <p'p'*> - d 2 <P ,2 > - (oa) 2 <u' 2 > 


p <P* > *• (Dfl) <U* >, 


For Gaussian distribution, f(k) * (“)^ 2 e /f(k)dk, we 

have p2 , p2 

,2 ,0, 3/2 12 1 . ,2 % 2 , 2 

<pt > „ (-*-) — - “ “j <p > « a <u' >, 

- p + P 

r l 2 P 


Quantities analogous to a square of the wave amplitude are 

ir 1 2 


A is A* -is' 

*** ' £/ S V s “ dE - rf i’ R ei e di: ' 


|?'| z Pj + p 2 . |GT - ijCPj + Pj), |p' 1^” p (Pj+p 2 ). 


I A I ^ ^ i _j - 5 

£ Mh? dk . - 2 t< k jp + < k 5 p 2 i. 

a k a 


|-.|2 . 1 


..2_ - 2 


A word is in order on the subject of units. The basic 

is is 

solutions A ■ (A /ak)R ,e n and P' - A R „e a have units of 
a a al a a a/ 


i 


J 


-1 7 ~ 

2-1 2 -9 

L T and L T *■, respectively, the same as the original variables, 

the velocity potential t and the pressure potential P. Thus the 

units of the amplitude A are I."t *" . Superposition of the basic 

solutions, their integrals over the wavenumber space, Eqs, (6), 

results in $ and P* becoming, respectively, the velocity potential 

density and pressure potential density, with units of L *T * and 

L ' l, T . Since f ■ |Ai", the units of the wavenumber distribution 

4 _4 

function f(k) are L T , and those of the probability density P 

are LT , However, the expectation values of any quantity will 

have units the same os those of the quantitv, thus the units of 

<u'*> are those of velocity square, and those of <p‘ > are pressure 

square. We may add that the velocity potential i- and the basic 

is 

solution 0 » (A /aklk ,e have units of action density as do 

a ot u 1 

the coefficients r, and f, . 

a n 

The root-mean-square amplitudes, due to a douhte integration 
over the wavenumber space under the square root sipn, should be 
interpreted as densities of wave properties, e,R. f /p r p r * * wave 
pressure density (units of pressure/volume), /u r u r * * induced 
velocity density, etc. 


< 


J 
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VI. i‘o*:i i.r.pns 


We have found that there exist situations In the acoustics 
of nonuniform media where classical treatment ho comes inadequate. 
To treat such situations quantum-like statistical methods were 
used to derive transport equations for the probability densities 
which determine wavenumber distribution functions. The distri- 
bution ''unctions permit one to evaluate ensemble avcrap.es as well 
as spectral distributions of physical quantities of interest. 

The present theory, it mav he argued (see, e.p., pp. 52-55 
of Xentzer, 1975b) , parallels that of quantum mechanics and 
oaeys the UNCEPTA 'NTY PRINCIPLE (existence of indeterminacy of 
phase or non-conservation of the number of waves), COMPLEMENTARITY 
PRINCIPLE (existence of pairs of conlugate variables, such as x 
and k, u> and t, each of which nay be better defined only at the 
expense of a correspond ine loss in the degree of definition of 
the other), and the CORRESPONDENCE PRINCIPLE (results of the 
theory reduce to the classical acoustics of homogeneous media). 

The nonlinear effects manifest themselves in the present 
theory in two ways. First, the hydrodynamic form of the transport 
equations contains nonlinear source terms and the nonunifonnity 
of the medium results in the coefficients of the transport equa- 
tions being functions of the dependent variables. Second, the 
bilinear terms introduce nonlinearity at the microscopic level 
through the wave interaction integrals. Thus each wave of finite 
amplitude affects amplitudes of all other waves. 

The fact that there exist two orthogonal (independent) wave 
solutions corresponding to a given wavenumber vector requires that 
a general solution or its integral representation be a aura of two 
Independent solutions (two modes of propagation), and that two 


I 


I 


I 


s 
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different sets of boundary conditions mnv be used for each mode. 
Thus it Is possible to impose a radiation condition at Infinity, 
and the coexistence of the two nodes at a ptven point may be 
viewed as a radiation condition at the microscopic level. 


u 


1! 


I 


-2M- 


te'VvmvK 


Pauli *3 lii,:cotiiiioal tlon of t he convective Have nutation 


Consider the potential equation. 


2 , 

+ 2u*'j'.“ + u* (u*’"'Vr; | - a * 0. 

St" 


This is equivalent to 

ff“r + u 0 ‘-rK^ + 


(nS.) 7 )* - « 


(A. 1) 


where the subscript ( ) Q denotes that u, is treated as a constant 

when operated upon hv p- 2 *. ** Trr* + u** 1 - 7 , and a 0 is kept constant when 

l' 8 t ot 

operated bv ? 0 . Equation (A.l) is the wave equation in the Fuler’s 
form for a fluid particle with an instantaneous velocitv u. 


In order to factor out the second order wave operator of 
Fq. (A. 11 into two linear operators, we write 

(i? - tT 2 )* » II, (A, 2) 

and require that 

6 <j> » 0 CA. 3) 

where E, Jl, and 0^ are linear differential operators. Then we 
need an operator fl., such that 

* (£ 2 - ft") * 0, 

The operators 0^ and 6 ^ must be real and must commute. If we choose 
Oj «* £ + H, and 0^ " E - H, 

than £ and i? must also commute. 


Following Pauli (see any book on the quantum theory) , we repre- 
sent £ and H by matrices and 4 by a vector. Taking 


E 


o 


D. 

4DTt 


a a 3 * 3 ^ 3 , 

H - a 3 ^U 
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wherc 5^,. . ,cr^ arc square matrices to be determined. We note that 
the operators E and H commute If u. and a. are treated as constants 

and if a, commutes with oL , eL and cL. Squaring ft and requiring 

'■2 ** 2 2 1 L 1 
that H “ , we have the following conditions 


«2 a 2 * -2 *■ 

a l ** a 2 ** a 3 “ 1 unit matrix), 

*1®1 ’ a 1^3 " * ot 3 Ct l* *2^3 “ *®3 Cl 2‘ 

Thus bt^, ct^, an ^ ^3 anticommute. Choosing - I, we aasure that 
commutes with S^, an ^ “ 3 * while the latter may be taken to 
be the Pauli spin matrices, e.g., 


0 1 


0 -i 


1 0 

1 u 

• *2 • 

i 0 

» S3 ~ 

0 -1 


Thus the solution of Eq, (A, 3) should be represented bv a two-component 
column vector, 



The system (A.3) is thus equivalent to 



+ 


a 



1 




2 


0 






0. 


(A.A) 


A* may be readily verified, as long as u. and a* are kept constant, 
or <t >2 may be eliminated by cross differentiation so that both 
components of ^ satisfy the convective wave equation (A. 1). With 
u» and a„ kept constant, Eqa. (A.A) are linear and we may consider 
$ to be a sum of two independent eigensolutiona. 



1 


0 

$ “ + * 2 * 4 *^ 

0 

+ < t >2 

1 


V*2 " °» 


and the system (A.A) decouples to yield 



IF 


f 


;; 


- 22 - 


J)«6. D*6- 36- 

— — +■ a* — " 0, and •*- — - a* — * 0, (A* 5) 

D,t 3z H,t 3* 

Bv permutations of the matrices n^ t we may single out 

any one of the three space coordinate* x, v, z, Thus a matrix 
permutation is equivalent to a rotation of the coordinate axes by 
90*. Since 36^/3* " n»V if n is a unit vector pointing in the 
direction of the preferred coordinate (z-axis in this case), then 
an arbitrary rotation of coordinate axes generalizes (A. 5) to 


D.6 

D.6 


* 4* a,n*?6^ “ [-i— *+• (u„+a»n> • V]6^ “ 0 


2 a.n*?4» 2 - + (u.-a.rj)»7]6 2 - 0 


CA.6) 


where n is arbitrary. The expresaiona ^ + (u,+a.n) *V in 

Eqs. (A. 6) are the total derivatives along the characteristic 
rays with n serving as a puraweter. Thus d/dt la an Eulerian 
derivative written for a nsuedo-particle (a wave packet) moving 
with a characteristic velocity (group velocity) u+an. 


Keeping n fixed implies that the solution vector 4 coincides 
with a different pair of eigensolutions at each point in time and 
space. If we single out one such a pair initially and would like 
to follow their development in time, we would have to determine 
how the unit vector changes in time. These two choices corres- 
pond to the Heisenberg's and Schrodinger's pictures of quantum 
mechanics. If necessary,, three equations for the components of 
the unit vector n(x,t) may be obtained as follows: 


n*n - 1, 


d(n»n ) „ 0 

dt 


and t|^ +(u-aH)«V][|^+(5+a5).V]6 a - I C§^) 2 -(a.?) 2 )6 a - 0. 

where a - 1,2. The unit vector n may be thus determined as a 
function of the solution of the wave equation. 


In order to derive the nonlinear convective wave equation, 

Eq, (A.l), all we need to postulate is that the dispersion relation 
for an arbitrary acoustic wave with a wavenumber vector k ■ kn at 
a point (x,t), where the local velocity and speed of sound are 
u(x,t) and a(x,t), is u » D in the inertial frame of reference 
which, at each point x « x, and time t « t. moves with different 
constant velocity u # + a 0 n equal to the local group velocity of 
the wave. Requiring that this postulate hold for an arbitrary 
wave (arbitrary n) anywhere (arbitrary x») and at all tiroes (arbi- 
trary t,), the dispersion relations, either w - u.*ic + a.n*£ * 0 

_ _ 2 2 2 

or (to - u, *k) - a,k “ 0, roust be viewed as a result of trans- 

forming ti) » 0 from an inertial frame instantaneously coinciding 
with the wuve packet and moving with the group velocity of the 
packet to an inertial frame common to all wave packets in all 
space-time. Thus u.(x,t) and a„(x,t) should be treated as para- 
metric constants when the dispersion relation is replaced by a 
differential operator, and should be allowed to become functions 

of x and t after the differential expression is expanded. Thus, 

- - 2 2 2 

e.g. , to the relation (u»-u,*k) - a.k m 0 there corresponds 

[i(~ + G.*V)] 2 $ + a 2 ? Z tJ> - 0. 

Expanding, we have 

2 

-(£- 4- + u.*?)* + + a 2 V 2 (f> 

at at D.tr 

2 

■ - + 2u.*|“-h u,* (u,*?)7<|>] + a 2 7 2 <J>- 0. 

Dropping the subscript ( ) , we obtain the nonlinear potential 
equation which, wh a written in terms of the operator D«/D«t 
takes the simple form of Eq. (A.l). 
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